GLOBAL WELL-POSEDNESS FOR THE 2 D QUASI-GEOSTROPHIC 
EQUATION IN A CRITICAL BESOV SPACE 
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Abstract. We show that the the 2 D quasi-geostrophic equation has global and unique 
strong solution, when the (large) data belongs in the critical, scale invariant space B^^ a D 

£2/(2a-l)_ 



1. Introduction 

In this paper, we are concerned with the mathematical properties of the Cauchy problem 
for the quasi-geostrophic equation in two spatial dimensions 

m 9 t + K(-A) a 9 + J(9)V9 = 

(l) 9(0, x) = 9°(x), 

where J (9) = (—R 2 9, R\9) and a £ [0, 1]. The physical meaning and the derivation of 
dU has been discussed extensively in the literature, we refer the interested reader to the 
classical book of Pedlosky, lfl2l . 

Depending on the value of the parameter a, one distinguishes between the subcritical 
case, a > 1/2, the critical case a = 1/2 and the supercritical case a < 1/2. It is known 
that the critical case a = 1/2 is especially relevant from a physical point of view, as it is 
direct analogue of the 3 D Navier-Stokes equations. On the other hand, considering the 
family of equations © with a £ [0, 1] allows us to understand better the influence of the 
diffusion on the evolution. 

An important scale invariance, associated with © is that 9 x (t, x) = \ 2a ~ 1 9(X 2a t, Ax) is 
a solution, if 9 is. It follows that the space H 2 ~ 2a (R 2 ) is critical for the problem at hand. 
A heuristic argument can be made that a well-posedness theory for initial data in H s , 
s < 2 — 2a should not hold. Thus, we concentrate our attention to the case s > 2 — 2a. 

The theory for existence of solutions and their uniqueness vary greatly, according to 
the criticality of the index a. For the crtical and supercritical case, the question has been 
studied in Q, 0, @, Q, lfT51 . |[Tr51l among others. The results are that when the data is 
large and belongs to H s , s > 2 — 2a, then one has at least a local solution, which may blow 
up after finite time. For small data in the critical space (or some Besov variant), Chae-Lee, 
1 1 ] and then J. Wu, lfl31 . ifTol have been able to show existence of global solutions. 

We would like to mention that the majority of these results have been subsequently 
refined to include Besov spaces of initial data with the same level of regularity and scaling 
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as the corresponding Sobolev spaces. Also, various uniqueness and blow-up criteria have 
been developed, see for example Section below. However, the fundamental question for 
existence of global, smooth solutions in the supercritical case remains open. We note that 
very recently, in the critical case a = 1/2, Kiselev, Nazarov and Volberg, [1 1J have shown 
the existence of global and smooth 1 solutions for any smooth (large) initial data. 

In the subcritical case, a > 1/2, which is of main concern for us, the quasi-geostrophic 
equation is better understood. Local and global well-posedness results, as well as LP decay 
estimates for the solution has been shown. 

To summarize the latest results, Constantin and Wu, [ 3 1 have shown global well-posedness 
for (the inhomogeneous version) of ©, whenever the data is in H s : s > 2 — 2a. For small 
data, there are plethora of results, which we will not review here, since we are primarily 
interested in the large data regime. On the other hand, time-decay estimates for 
have been shown in Q and 0, see Section[2]below for further details. Finally, we mention 
a local well-posedness result for large data in H 2 ~ 2a fl L 2 , due to Ning Ju, J9|. Note that 
the space H 2 ~ 2a is not scale invariant (due to the L 2 part of it) and thus, such solutions 
cannot be rescaled to global ones. 

In this work, we show that the quasi-geostrophic equation is globally well-posed in the 
critical space B^lx* ^ L 2 ^ 2a ~ l \ that is whenever the data 9° belongs to the space, there is 
a global and unique 2 solution in the same space. 

Theorem 1. Let a G (1/2, 1). Then for any initial data 9° G B 2 2 ~ 2a {K 2 ) n L 2 ^ 2a ^(R 2 ), 
the quasi-geostrophic equation © has a global solution 

9 G L°°([0, oo); B^ a (R 2 ) n L 2 ^ 2a_1 ^(R 2 )) 

Moreover, the solution satisfies the a priori estimate 

(2) ||^(^)||^2-2 QlnL 2/(2a-l) < C Kja (||6 |0 |L2-2a ni 2/(2 a -l) + 1 1 9° 1 1 2 /(2a-l) ) 5 

for allt>0 and M(a) = max(2, l/(2a — 1)). In particular, the norms remain bounded 
for < t < oo. 

In addition, if9 G L 2 (R 2 ), then 9 G L 2 ((0, oo), H a (R 2 )), in fact 

(3) II^IIl 2 ,//"(r 2 )) — II^°IIl 2 (r 2 )- 

For a fixed T > 0, the solution is unique class of weak solutions on [0,T] satisfying 
9 G L°°([0, T], L 2 (R 2 )) fl L°°([0, T], L 2 ^ 2a ^) fl L 2 ((0, T), H a (R 2 )). 

Several remarks are in order. 

(1) Note that global solutions exist and are unique in the space i?2^ Q (R 2 )nL 2/(2a_1) (R 2 ), 
when the data is in the same scale invariant space. Note that such space properly 
contains H 2 ~ 2a (R 2 ). In other words, taking data in H 2 ~ 2a (R 2 ) guarantees the ex- 
istence of global solution, but by © we only know that the slightly smaller norm 

o2-2c n r2/(2a-i) stays bounded. 

2,oo 

lr The smoothness assumption in II II is essentially at the level of H 2 (R 2 ), while the the critical case, the 
critical Sobolev space is iJ 1 (R 2 ). 

2 For the uniqueness one has to assume in addition 6>° 6 L 2 (R 2 ) 
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(2) It is an interesting question, whether Theorem [T] and more precisely © holds in 
in the case of the Sobolev space H 2 ~ 2a or even for some Besov space in the form 
B 2 ~ 2a for some r < oo. We note that the main difficulty is the estimate © for 
smooth solutions, since once © is established, one easily deduce the global exis- 
tence and uniqueness by standard arguments. 

(3) The results in Theorem [l] apply as stated for the case of the quasi-geostrophic 
equation on the real line, but one can recover the exact same statement, if one 
considers on T 2 . We omit the details, as they amount to a minor modification 
of the proof presented below. 

Acknowledgement: I would like to thank Ning Ju for several stimulating discussions on 
the topic. 

2. Preliminaries 

2.1. The 2 D quasigeostrophic equation - existence and maximum principles. We start 
this section by the Resnick's theorem, lfT3l for existence of weak solutions. That is when- 
ever 9° G L 2 (R 2 ) and for any T > 0, there exists a function 

9 G L°°([0, T], L 2 (R 2 )) n L 2 [[0, T], H a (K 2 )), so that for any test function <p 

[ 9(TMT)- [ T [ 6{J{6)V<p) + K f [ ((-A)^)((-A)^)=/^(0,4 
Jr 2 Jo Jr 2 Jo Jr 2 J 

In his dissertation, JT3l . Resnick also established the maximum principle for L p norms, 

that is for smooth solutions of © and 1 < p < oo, one has 

(4) ||0(f) || 

This was later generalized by Constanin-Wu, Q, for the case p = 2 and by Cordoba- 
Cordoba, [ 5 1 in the case p = 2 n and N. Ju, (7 ] for all p > 2 to actually imply a power rate of 
decay for ||^(t)|| iP ( R 2) and an exponential rate of decay, when one considers the equation 
on the torus T 2 . In the sequel, we use primarily ©, but is nevertheless interesting 
question to determine the optimal rates of decay for these norms. Note that Constantin and 
Wu have shown in Q, that the optimal rate for \\9(t) \\ L 2m^ is < t >~ 1 / 2a . Ning Ju has 
proved in [7|, that 3 ||0(t) || LP(R2) < C(||0°|| LP )(1 + t)-^)^. 

2.2. The uniqueness theorem of Constantin- Wu. Recall the uniqueness theorem of 
Constantin-Wu (Theorem 2.2, in 0) 

Theorem 2. (Constantin-Wu) Assume that a G (1/2, 1] and p, q satisfy p > 1, q > 1 and 
1/p + a/q = a — 1/2. Then for every T > 0, there is at most one weak solution of © in 
[0, T], satisfying 

9 G L°°([0,T],L 2 (R 2 )) nL 2 [[0,T],H a (R 2 )) n L 9 ([0, T], L P (R 2 )). 

In particular, one can take q — oo, 1/p — a— 1/2 to obtain uniqueness for weak solutions 
satisfying 9 G L°°([0, T], i7(R 2 )). 

3 For example, (p — 2) /2pa — > as p — > 2, whereas the optimal rate is (2a) _1 , as shown by Constantin 
and Wu. On the other hand, we must note that the rate of L p decay obtained by Ning Ju holds under the 
assumption that 9 a 6 L 2 (R 2 ), while Constantin-Wu assume that 9q G L 1 (R 2 ). 
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2.3. Some Fourier Analysis. Define the Fourier transform by 

.AO = J f(x)e- ix <dx 

and its inverse by 

R n 

For a positive, smooth and even function % : R 1 ^ R 1 , supported in {£ : |£| < 2} and 
x(0 = 1 for all |£| < 1. Define </?(£) = x(0 ~~ x(2£)' which is supported in the annulus 
1/2 < ICI < 2. Clearly £ fee2 = 1** all ^ 0. 

The k th Littlewood-Paley projection is P k f(g) = ^( 2_fe 0/(0- Similarly P <k = £ P, 

/<fc 

given by the multiplier x(2~ fc £) etc. Note that the kernels of P&, P <k are uniformly in- 
tegrable and thus P k , P <k : L p — > L p for 1 < p < oo and ||Pfc|| L p_^ LP < C||x|| L i- m 
particular, the bounds are independent of k. 

The kernels of P fc are smooth and real- valued 4 and P k commutes with differential oper- 
ators. We will frequently use the notation ip k (x) instead of P k t/j, when this will not create 
confusion. 

It is convenient to define the (homogeneous and inhomogeneous) Sobolev norms in 
terms of the Littlewood-Paley operators. Namely for any s > 0, define for every Schwartz 
function ip th norms 



£ 2 2ks \\A 



1/2 




and the corresponding spaces are then obtained as the closure of the set of all Schwartz 
functions in these norms. Clearly H s = L 2 n H s . 

Introduce the operator A acting via A?/>(£) := |£|V>(£). Clearly, by the uniform bound- 
edness of P k in the scale of L p spaces, ||A s *0 fe || L j, ~ 2 fcs ||^ fc || LJ( . 

Next, we introduce some basic facts from the theory of the paraproducts, which will be 
useful for us, when estimating the contribution of the nonlinearity. 
Write for any two Schwartz functions /, g and any integer k, 

Pk(fg) = Pk(52fh9h) = Pk( £ fhQh) + Pk{ £ f h g h ) 

h,h h,h-\h~h\<3 h,h:\h-h\>3 

But 

p k( £ fi^h) = p k{ £ hsh) 

h,h-\h-h\<3 h,h:\h-h\<3,min(h,h)>k-3 



4 Thus for a real valued function ip, Pkip is real-valued as well. 
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since by the properties of the convolution 2 h+1 + 2 l ' 2+1 must be at least 2 fe_1 and 

(l,^2:|2l— ^2 1 >3 h,h'-\h— ^ 2 1 > 3, | max(Zi,Z2)— k\<3 

since otherwise suppf^g^ C {£ : |£| ~ 2 max ( il >' 2 )}, which would be away from the set 
~ 2 fc } and thus P k (f h g h ) = 0. 
All in all, 



(5) 



Pk(fg) = P k ( E PifPi-3<-<i+3g)+ 

l=k-3 

3 3 

Pk{ Pk+jfP<k+j-3g) + Pk( J2 Pk+jgP<k+j-3f)- 

i=-3 j=-3 



We will refer to the first term as "high-high interaction" term, while the second and the third 
terms represent the "high-low interaction" term. We have the following lemma, which is 
an application of the representation formula ©. 

Lemma 1. For every < s < 1, 2 < p, q < oo : 1/p + 1/q — 1/2, there is the estimate 
| / P k M-JW ■ W]dx\ < C s 2 k ^\\ip k \\ LP { 2- s(i - fe) ||A^z| 

J l>k-3 

for some absolute constant C. 
Proof. Integration by parts and div(J(6)) = yield 

PkMJW ■ Vipjdx = - [ WiP k P k • [J(iP)iP]dx 



\l 2 



Li- 



At this point, by the boundedness of the Riesz transform on LP, we treat J(ip) as Tip, where 
T : U — > U for all 1 < r < oo and ignore the vector structure. By Holder's inequality, 



ViP k P k [T(iP) ■ iP}dx\ < 2 fc ||^|| LP ||P fc [T(^]|| 



LP 1 ' 



By ©, 

\\P k [T(iP)iP]\\ LP , < 



Y PiTipPi- 3 <.<i +3 ip 



l=k-3 



LP 



Y Pk+j(Tip)P <k+ j-3ip 



i=-3 



+ 



LP 



+ 



LP 



l=k-3 



Y Pk+j(^P)P<k+j-3Tlp 
3 

Y \\P k+ ML4P<k +j -3^\\ Lq <C(Y, 



< Y IWIUI^3<.<1+3^IIl« + 



>l\\V> 



Li- 



i=-3 



l>k-3 



The Lemma follows by the observation HV'zll^ ~ 2 ||A S ^/|| L2 and by reshufling the 2 . 



□ 
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3. Proof of Theorem [T] 

The main step of the proof of Theorem [T] is the energy estimate ©. 

We start with the assumption that we are given a smooth solution 9{t, x), corresponding 
to an initial data 9° up to time T and we will prove © based on it. Assume © for a moment 
for such smooth solutions. We will show that the global existence and uniqueness follows 
in a standard way from an approximation argument and the Constantin-Wu uniqueness 
result, Theorem^ 

Indeed, for a given initial data 9°, take an approximating sequence in B\ n L 2 ' ( 2a-1 ), 
{d®} of smooth functions (say in the Schwartz class S). By the Constantin-Wu existence 
result for data in H s : s > (2 — 2a), we have global and smooth solutions 9i(t). In 
addition, they will satisfy the energy estimate ©. Moreover, by the LP maximum principle, 
< ||0j(O)|| £9 for all 1 < q < oo, in particular for q = 2,q = 2/ (2a - 1). 

Taking weak limits will produce a weak solution 9{t) of ©, corresponding to initial 
data 9°, so that it satisfies the energy estimate © and ||0|| LOO £2/(2a-i) < ||6> || L 2/(2 Q -i). This 
shows the existence of a weak solution with the required smoothness of the initial data. 

For the uniqueness part, we should require in addition that 9° E L 2 (R 2 ). Then, we 
will show || 1| < ||#°|Il2, which allows us to apply the Constantin-Wu uniqueness 
result (Theorem That is, 9 is the unique solution in the class L°°([0, T), L 2 (R 2 )) n 
L 2 [[0, T], £P(R 2 )) n L°°{[0, T], L 2 ^ 2a ~^ (R 2 )). Thus, it remains to prove © for smooth 
solutions and ©. Since, © is relatively easy, we start with ©. 

3.1. Proof of the energy estimate ©. Let s = 2 — 2a. Take a Littlewood-Paley operator 
on both sides of © 

d t 9 k + n(-A) a 9 k + P k (J(9)V9) = 0. 
Taking a dot product with 9 k (which is real-valued!) yields 

d t \\9 k \\ 2 L2 + 2/ t ||(-A) Q / 2 ^||^ 2 + 2 j P k 9 k J(9)V9 = 0. 

By the properties of the Littlewood-Paley operators, ||(— A) a / 2 0fc||^ 2 ~ 2 2afc ||6'fc|| 2 : 2- For 
the integral term, use Lemma[T]with 1/p = 1/2 — sq/2, 1/q = s /2. We have 

| / P k 9 k J{9)V9dx\ < C2 k{1 -^\\9 k (t)\\ Lq { 2" so( ^ ) ||A^|Il 2 )II^(0IIl, < 

l>k-3 

<C2 k ^\\9 k (t)\\ LqS np\\A%\\ L 49(m LP - 

i 

By the LP maximum principle, @, we have ||#(£)|| LP < Substituting everything 

in the equation allows us to conclude 

(6) dtPk\?v + cK2 2te ||0 fc ||J a < C72 fc ( 1 -" > )||0 o || LP ||0 fc (t)|| £ ,sup||A»°0 l || La 

At this point, the argument splits in two cases with a treshold value of a = 3/4. As 
expected, the case 3/4 < a < 1 proves out to be slightly simpler, so we start with it. 
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3.1.1. The case 3/4 < a < 1. The significance of the restriction a > 3/4 is in the fact 
that s = 2 - 2a e (0, 1/2]. Therefore 1/q = s /2 < 1/2 - s /2 = 1/p, implying p < q. 
Thus, by the Sobolev embedding 5 , the boundedness of P k on L p and the L p maximum 
principle imply 

\Mt)\\ Lq < 2 2k ^- 1 ^\\9 k (t)\\ LP < 2 fc ( 1 - 2 ^||^(t)|| iP < 2 fe ( 1 - 2s °)||0°|| iP 
By ©, we infer 

(7) d t \\9 k \\ 2 L2 + cK2 2ka \\9 k \\ 2 L2 <C2 fc ( 2 - 3a °)||0 o ||^sup||A flo ^(*)llL 2 

It is a standard step now to make use of the Gronwal's inequality, namely rewrite as 

and estimate after integration 

(8) \\e k (t)\\ 2 L2 < C^-^-^pm sup sup \\A s %(z)\\ L2 + \\9° k \\ 2 L2 e-™ 22kat . 

0<z<t I 

Note that in the formula above C k ~ 1/k and 2 — 3s — 2a = —2s . 
Introduce the functional 

J(t) = sup su P 2 feS0 ||^(2)|| L2 . 

0<z<t k 

Clearly, one may deduce from © that 

J 2 (t)<J 2 (0) + C K J(t)\\9°\\ 2 LP , 

hence 

j(t)<2j(o) + c K \\e°\f LP , 



which is 

(9) su P 2 fc ( 2 - 2Q )||^(t)|| L2 < 2sup2 fc ( 2 ~ 2Q )||^|| L2 +C K \\9°\\ 2 LP . 

k k 

This is the a priori estimate of the solution 9, © for the case a G [3/4, l). As we have 
observed in the beginning of the section, it follows that the 2 D quasi-geostrophic equation 
([l]) has global solution with (potentially large) data in the scale invariant space B 2 ^ a (R 2 )n 
L 2 /(2«-i)(R2). 

3.1.2. The case 1/2 < a < 3/4. In this case, it is clear that s = 2 - 2a e (1/2, 1), 
whence 2 < q = (1 — a)~ l < p = (a — 1/2) -1 . By Gagliardo-Nirenberg's, 

||0fc||z<! < c||a 2 2Q #fc||^ 2 ||A -a ^ILp 7 ' 

with 7 = §5|f e (0,1) and a = (2 ~ 2 ^i 3 ~ 4a) . Thus, by \\A- a 9 k \\ LP ~ 2- ak \\9 k \\ LP , whence 
it follows that 

\\o k \\ Lq < c2- fc ( 3 - 4Q )sup||A so ^n^ii4Wiii; 7 . 

i 



5 or more appropriately the Bernstein inequality 
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Substituting this in © yields 

do) d t \\e k (t)f L2 + c K 2 2ka \\e k (t)f L2 <2 fc ( 1 - s ''- 3+4Q )sup||A s ^niit 7 II^WIlL,||^°llL7- 

Using the maximum principle ||0fc(t)||£i> < ||^°||ip, this reduces to 

d t \\9 k (t)f L2 + CK2 2ka \\9 k (t)\\l 2 < 2 fc ( 1 - s «- 3+4Q ) sup ||a so ^||^||^ || 2 .; 7 . 



By the Gronwall's inequality, we deduce 



(11) \\0 k (t)\\ z L2 < ^ > 1 + C K 2^ S0 sup sup ||A s °0 z (z)||£ 7 

0<z<t I 



By using the same energy functional J(t) defined above, we conclude that 

J 2 (t)< j 2 (o) + c k [j(t)] 1 ^\\e°\\ 2 LP \ 

Since 1 + 7 < 2, by Young's inequality 

^(t)<J 2 (o) + ^ + a,j^iii 4 ; 27)/(1 - 7) . 

whence 

J(t)<2J(0) + C K>J \\f)' 

which is 

(12) su P 2 fe ( 2 - 2Q )||^(t)|| L2 <sup2 fc ( 2 - 2a )||^|L 2 + a, 7 ||^ 



0||(2- 7 )/(l-7) 
LP 



0|| (2- 7 )/(l-7) 



L 2 ' -K,7ir 1 1 lp 

K K 

Again, this implies © with M(a) = 1/ (2a — 1) and the problem ^ has global solution 
in ^^(R 2 ) n £ 2/(2a ~ 1} (R 2 )> when the initial data is taken in the same space. 

3.2. 9 G L oo ([0,cx)),L 2 (R 2 ))nL 2 ((0,oo), J ff a (R 2 )). Both of these estimates are classical 
for smooth solutions, but we sketch their proofs for completeness. 

In fact, 6 G L°°([0, 00), L 2 (R 2 )) follows from the maximum principle ©. For the 
second estimate, we multiply the equation by 9 and integrate in x. We get 



dt\\e(t)\\l* + \\^ a o(t)\\h = -] 0[J(0)V6}dx = 

Time integartion now yields 

riiA^ii^^^ii^ii^-ii^iii^ii^ii 2 ,, 

Jo 



whence 6 G L 2 ((0, 00), H a (R 2 )). 
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